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Abstract 

The three different hehcity states of bosons, produced in the reaction e"'"e~ 
W"'"W~ £z/qq are studied using leptonic and hadronic W decays at ^/s=183 GeV 
and 189 GeV. The W polarisation is also measured as a function of the scattering 
angle between the W~ and the direction of the e~ beam. The analysis demonstrates 
that W bosons are produced with all three helicities, the longitudinal and the two 
transverse states. Combining the results from the two center-of-mass energies and 
with leptonic and hadronic W decays, the fraction of longitudinally polarised 
bosons is measured to be 0.261 ± 0.051(stat.) ± 0.016(syst.) in agreement with 
the expectation from the Standard Model. 
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Introduction 



Previous measurements of W"'"W~ production at LEP have concentrated on measurements of 
the W mass, the W branching ratios, the differential and total cross sections and the anomalous 
couplings PHHI- These measurements show, using the differential cross sections with respect to 
the W production and decay angles, good agreement with theoretical calculations within the 
Standard Model [HE]- This good agreement with the Standard Model indicates indirectly that 
W bosons with all three helicities are produced in the reaction e"'"e~ —>■ W+W~. 

The primary goal of the measurement described in this paper is a quantitative and model 
independent analysis of all three W helicity states and in particular, the direct observation of 
longitudinally polarised W bosons. Measurements of longitudinally polarised W bosons have 
previously been reported in the reaction e+e" W"''W~ and in top decays jS] • 

At center-of-mass energies close to 190 GeV and within the Standard Model, one expects 
that about one quarter of all W bosons should be longitudinally polarised [7j. Furthermore, 
the production of W bosons with different helicities depends strongly on the W~ scattering 
angle 6'w- with respect to the e~ beam direction. For example one expects for larger 
than 90 degrees that almost 40% of the events contain at least one longitudinally polarised W 
boson. In contrast, for 6'w- between 20 and 70 degrees, the cross section is dominated by the 
neutrino-exchange diagram and the W"'"W~ should be produced dominantly with transverse 
polarisation. The fractions of the helicity states should thus also be measured as a function 
of ^w- • 

The measurement is performed with the L3 detector at LEP, using data samples of 55.5 pb~^ 
and 176.4 pb~^ collected at average center-of-mass energies of 183 GeV and 189 GeV, respec- 
tively. A detailed description of the L3 detector and its performance is given in reference [H]. 
The L3 detector response for W+W" events from the KORALW [9J and the EEWW [lU] Monte 
Carlo programs is simulated with the GEANT-based L3 detector simulation program jllj . 

Analysis strategy 

The different W helicity states result in different angular distributions of the W decay prod- 
ucts. The decay angle 6* in the W rest frame between the left-handed negatively charged 
lepton and the W~ has a (1 ± cos^*)^ distribution for a W~ with helicity =f1. The right- 
handed positively charged lepton has a (1 ± cos 6'*)^ distribution for a W"*" with helicity ±1. 
Longitudinally polarised W bosons (helicity 0) result in a symmetric distribution of the decay 
products, proportional to sin^O*. To simplify the description of the helicity fractions, we refer 
in the following text only to the fractions /_, /o and /+ of the W~ helicities, which includes 
the corresponding W"*" states with /+, /o and /_, respectively. 

In order to study the W polarisation, we use events of the type e"'"e~ W^W^ £^z/qq 
with being either e^ or fi^. The neutrino four-momentum vector is reconstructed from the 
total missing momentum vector of the event. These event samples are essentially background 
free and allow a measurement with good accuracy of the momentum vector, the W charge 
and the decay angle 6* in the W rest frame. 

In contrast to leptonic W decays, where the decay angle 9^ of the £^ is well defined, the 
corresponding 6* for quarks in W decays has to be calculated from the hadronic decay prod- 
ucts. To approximate the quark decay angle in the W rest frame, we proceed in the following 
way. First, all particles besides the charged lepton and the missing neutrino in the event are 
associated with the hadronic decay of the W. We then calculate their associated four-vectors 
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in the rest frame of the W and determine the corresponding thrust axis in this rest frame. The 
angle 0^-^^^^^ of this thrust axis with respect to the W momentum vector in the laboratory frame 
is used to describe the quark decay angle 9* in the W rest frame. 

After correcting for efficiencies, the contributions from different W polarisation states are 
obtained from a fit to the cos^* distributions. For the leptonic W decays the fractions /_, /+ 
and /o of the three W helicity states are obtained from: 



fj- (l + cosn' + /4 (l-cosn^ + Zo^in^r. (1) 



NdcosO* ' 8 ' ' '^S ' /I 

For hadronic W decays, without quark charge identification, one measures only the absolute 
value of the W hadronic decay angle | cos Q* \ . However, this distribution can still be used to 
measure the fractions for the sum of the two transverse helicity states f± = f- + /+ and /o 
using: 



/J(l + cos2r) + /o^sin2r. (2) 



Nd\cose*\ '^4' 

The predictions for the compositions of W helicity states as a function of the W~ scattering 
angle 6^^- , following the formalism of Hagiwara et al. [7] and its implementation in the KO- 
RALW Monte Carlo program P, are used as the Standard Model prediction for our analysis. 
The helicity composition of the total W sample is extracted from a fit to the distribution of the 
simulated decay angles. From a fit to a KORALW Monte Carlo event sample at y^=189 GeV, 
with a size 100 times larger than the data sample, the Standard Model predictions for inclusive 
W helicity fractions /_, /+ and /o are obtained to be 56.3%, 18.0% and 25.7%, respectively. 
The statistical errors are smaller than 0.5%. 

Within the statistical errors, the same fractions are found from a WW event sample gener- 
ated with the EEWW Monte Carlo program PO] which uses the zero total W width approxima- 
tion and assigns the W helicities on an event by event basis. The W helicity fractions obtained 
from the fit to the decay angle distributions agree, within statistical errors smaller than 0.9%, 
with the generated W helicity fractions. This shows that the Born level formulas (1) and (2) 
are applicable after radiative corrections. 



Selection of W+W — > e(/x)i>'qq events 

The selection of W+W~ — > e(/i)i/qq events is similar to the selections described in our previous 
publications on WW final states P]. However, in order to assure well measured W production 
and decay angles, more restrictive criteria are used. Charged leptons are identified using their 
characteristic signatures. Electrons are identified as isolated energy depositions in the electro- 
magnetic calorimeter with electromagnetic shower shape which are matched in azimuth to a 
track reconstructed in the central tracking chamber. The energy and direction of electrons are 
measured using the electromagnetic calorimeter, while the charge is obtained from the associ- 
ated track. Muons are identified and measured with tracks reconstructed in the muon chambers 
which point back to the interaction vertex. All other energy depositions in the calorimeters are 
assumed to originate from the hadronic W decay. The neutrino momentum vector is set equal 
to the total missing momentum vector of the event. In addition the following criteria are used 
for the selection of W^W~ e(/i)z/qq events: 
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• The reconstructed momentum should be greater than 20 GeV for electrons and 15 GeV 
for muons. 

• The momentum of the neutrino should be greater than 10 GeV and its polar angle, 9^, 
has to satisfy | cos6',^| < 0.95. 

• The invariant mass of the Cv system should be greater than 60 GeV. 

• The invariant mass of the hadronic system should be between 50 and 110 GeV. 

Using these criteria, 81 and 288 events of the type W^W^ ez/qq are selected at center- 
of-mass energies of 183 GeV and 189 GeV, respectively. The corresponding event numbers for 
/xz/qq are 67 and 262 events. Adding the electron and muon event samples together, we find 68 
and 280 i~ events and 80 and 270 events, respectively, in the 183 GeV and 189 GeV data 
samples. These samples have a purity of 96%, where the background from W"'"W~ — > ri^qq 
with leptonic r decays and the background from e"'"e~ — *■ hadrons contribute each about 2%. 

The measured cos 6'-w- distribution is found to be in good agreement with the MC expec- 
tations, as shown in Figure ^ for events with electrons and with muons for the 189 GeV data 
sample. About 5% of the accepted events with electron candidates have a wrongly assigned 
charge. Charge confusion is insignificant for events with muons. The charge confusion depends 
on the reconstructed W~ scattering angle and is largest for W bosons with small scattering an- 
gle with respect to the beam direction. This results in a small misassignment between W bosons 
with helicity +1 and -1 but has negligible effects for the fraction of longitudinally polarised W 
bosons, which is essentially independent of the charge assignment. 

Analysis of the W helicity states 

After subtracting the backgrounds from the data, the fractions of the W helicity states are mea- 
sured from the distributions dN /d cos 9^± and dN / d\ cos 9^^^^^^^] for the leptonic and hadronic 
W decay angle and as a function of the scattering angle 6'w- ■ 

To extract the W helicity fractions, the observed distributions are corrected for the selection 
efficiencies which are obtained as a function of cos 6'*. To take into account possible deviations 
between the helicity fractions in the data and Monte Carlo as a function of cos ^w- ) the data 
are corrected differentially using 9 bins of the cos 9^^- scattering angle. For each cos ^w- bin, 
the efficiency is obtained as a function of cos^^* using the ratio of the reconstructed and the 
generated cos 6** distributions for the leptonic and hadronic W decays. The measured cos6'* 
distributions for the corresponding cos ^w- bins in the data are corrected and combined. 

The efficiency corrections are obtained from large samples of fully simulated KORALW 
Monte Carlo events. Using these Monte Carlo events we have studied the accuracy with which 
we reconstruct the 9* decay angles. The study shows that 9* is reconstructed with a standard 
deviation of 9.2 degrees and a small shift of —3.2 degrees for the leptonic W decays. For 
hadronic W decays one finds that 9* is reconstructed with a standard deviation of 12.0 degrees 
and a shift of +3.3 degrees. 

The bias and sensitivity loss due to the efficiency corrections and the 9* resolution has been 
determined with fully simulated and reconstructed Monte Carlo events where the generated W 
helicity fractions have been varied over a large range. This was done both with the EEWW 
Monte Carlo program, where the generated W helicities are known on an event by event basis 
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and with the KORALW Monte Carlo using a weighting method to assign the W hehcities on a 
statistical basis, ignoring W spin correlations. 

Averaging both Monte Carlo estimates one finds that leptonic W decays with 100% helicity 
— 1 states would be measured to consist of 94% of helicity —1 and 6% helicity states while 
a W sample with 100% helicity +1 would be reconstructed to consist of more than 99% of 
helicity +1 states. Similar numbers are found if one starts with 100% helicity 0, which would 
be measured with 92% helicity , 3% helicity —1 and 5% helicity +1. The corresponding 
numbers for hadronic W decays are that 94% of W bosons with helicity ±1 and 85% of W 
bosons with helicity are correctly reconstructed. The study has been repeated as a function 
of cos ^w- and within the statistical errors the results are the same as the ones from the total 
W sample. To obtain a correction function for the bias and the efficiency loss, the fraction /q 
has been varied between and 100%. A linear relation between the generated and the fitted 
W helicity fractions is found. 

Results and systematics 

These efficiency corrected cos 9* distributions are used to extract the helicity fractions. The 
results of the binned fits to these distributions for leptonic and hadronic W decays from the 
ys=189 GeV data are shown in Figure El No constraint on the total cross section is applied 
and one finds that the data are well described only if all three W helicity states are used in the 
fit. Fits which include only —1 and ±1 helicities, as also shown in Figure 2, fail to describe 
the data. For leptonic W decays one finds that the increases from 7.1 for seven degrees of 
freedom if all three W helicity states are included to 17.8 for eight degrees of freedom if only 
helicity —1 and +1 are used to describe the data. For hadronic W decays the increases 
from 9.8 for eight degrees of freedom if all three W helicity states are included to 26.4 for nine 
degrees of freedom if only helicity ±1 are used to describe the data. 

The fraction of longitudinally polarised W bosons in the y^=189 GeV data is measured to 
be 0.220±0.077 for the leptonic decays and 0.285±0.084 for hadronic decays. The fractions for 
the different W helicity states, together with the Standard Model Monte Carlo expectations, 
are given in Table 1 for the i/s=189 GeV and A/i=183 GeV data. The observed fractions of 
longitudinally polarised W bosons measured with leptonic and hadronic W decays agree with 
each other and with the Standard Model expectation of 0.26 and differ from zero by several 
standard deviations. 

Systematic studies have been performed to verify the stability of the fit results with respect 
to the fraction of longitudinally polarised W bosons. We have investigated (1) uncertainties due 
to backgrounds, (2) efficiencies and selection criteria, (3) the hadron energy response functions 
of the electromagnetic and hadronic calorimeters, (4) the difference between the differential and 
overall efficiency corrections and (5) a method where the fraction /o has been obtained directly 
from a fit to the measured cos 9* distributions using the Monte Carlo shape from the different 
W helicity states after the reconstruction. 

The analysis has been repeated assuming large relative background uncertainies of ± 50% 
from either the hadronic background or from misidentified W— > tv decays. Using these mod- 
ifications the measured fractions of longitudinally polarised W bosons is found to vary by at 
most 0.012 for leptonic W decays and by 0.004 for the hadronic W decays. The hadron energy 
measurement is obtained from a combination of the energy deposited in the electromagnetic 
and hadron calorimeter multiplied by calibration constants which take the average calorimeter 
e^ /hadron response function into account. These calibration constants have been varied over 
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a wide range while demanding that the average of the reconstructed masses for leptonic and 
hadronic W decays agree within better than ± 3 GeV with an average W mass of 80.4 GeV. 
Since the neutrino momentum vector is reconstructed from the observed missing momentum 
vector, correlations exist between the reconstructed decay angles in the hadronic W system and 
the corresponding leptonic W system. For example, a particular choice of the energy calibration 
constants reduces the fraction of longitudinally polarised W bosons by 0.024 as seen with the 
leptonic W decays but increases the corresponding fraction for the hadronic decays by 0.015. 

Similar variations in the fraction of longitudinally polarised W bosons have been seen with 
the other systematic studies, as summarised in Table 2. Assuming that the variations given in 
Table 2 are all due to systematics and adding them in quadrature, a systematic error of ±0.034, 
±0.024 and ±0.016 is assigned to the fraction of longitudinally polarised W bosons measured 
with leptonic, hadronic decays and for the combined measurement, respectively. 

Combining the results from the ^/s = 183 GeV and 189 GeV, ignoring the slight energy 
dependence of the W helicity fractions expected from the Standard Model, the fraction of 
longitudinally polarised W bosons is measured to be 

/o = 0.261 ± 0.051(stat.) ± 0M6{syst.) 

and agrees with expectation from the Standard Model of 0.26. 

As mentioned in the introduction, it is interesting to measure the W helicity fractions as 
a function of the W~ scattering angle 6'w-- Thus the fits are repeated for different ranges of 
cos 6'-w-- The cos 6'-w- ranges are selected such that the contributions from the transversely 
polarised W bosons should be either suppressed or enhanced as shown in Figures El and El 

To obtain quantitative numbers for the W helicity fractions as a function of cos6'w- the 
data from the two different center-of-mass energies are combined and the helicity fractions 
are measured for three bins of cos^w The bins are chosen such that large variations of the 
different helicity fractions are expected [2] yet keeping a sufficient statistical significance. The 
results, given in Table 3, agree with the Standard Model expectations and demonstrate that the 
fraction of W bosons with helicity —1 depends on the W scattering angle as shown in Figure El 

In summary, all three W boson helicity states, the two transverse as well as the longitu- 
dinal ones are observed with fractions in agreement with Standard Model expectations. The 
production of longitudinally polarised W bosons is thus directly observed with a significance of 
five standard deviations. 
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^/i = 189 GeV 




Hc4icit}- AA^ ^ (ly 


Helicity AA"^ liadroiis 




-1 


+1 





±1 





Data 


0.568 ± 0.071 


0.212 ± 0.046 


0.220 ± 0.077 


0.716 ± 0.086 


0.285 ± 0.084 


MC 


0.56 


0.18 


0.26 


0.74 


0.26 


= 183 GeV 


Data 


0.56 ± 0.14 


0.10 ± 0.08 


0.34 ± 0.15 


0.75 ± 0.17 


0.25 ± 0.17 


MC 


0.53 


0.20 


0.27 


0.73 


0.27 



Table 1: Measured W helicity fractions for the leptonic and hadronic W decays for 
the \/s = 189 GeV and ^/s = 183 GeV data samples. The corresponding helicity 
fractions in the Standard Model as implemented in the KORALW Monte Carlo 
program where the statistical errors are negligible in comparison with the data are 
also given. 
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Fraction of longitudinally polarised W bosons 






W— > hadrons 


average 


standard method 


0.220 ± 0.077 


0.285 ± 0.084 


0.252 ± 0.057 


background corrections 
efficiency uncertainty 
calorimeter calibration (liadrons) 
integrated efficiency correction 
analysis method 


0.209-0.232 

0.214 
0.197-0. 215 

0.233 

0.237 


0.282-0.286 

0.279 
0.282-0.300 

0.268 

0.279 


0.241-0.258 

0.247 
0.244-0.254 

0.250 

0.258 



Table 2: Measurements of the fraction of longitudinally polarised W bosons for 
leptonic and hadronic W decays from the y/s = 189 GeV data sample investigating 
various sources of systematics. 
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Vi = 183 + 189 GeV data 




H(4icity ^ ( 




H(4icity AA^ 


—f luulroiis 


\AJ\j L/\Y~ 


-1 


+1 

1 J- 


n 


+1 

_l L 


n 






0.45 ± 0.22 




n 07 _|_ n oo 


n 1 Q _|_ n 00 


-0.4 - 0.3 


0.40 ± 0.09 


0.23 ± 0.08 


0.37 ± 0.12 


0.94 ± 0.16 


0.06 ± 0.15 


0.3- 1.0 


0.66 ± 0.08 


0.08 ± 0.04 


0.26 ± 0.08 


0.75 ± 0.11 


0.23 ± 0.10 






= 183 + 189 GeV KORALW MC 






Helicity W ^ 


Helicity W 


hadrons 




-1 


+1 





±1 


n 




0.13 


0.45 


0.42 


0.58 


AO 


-0.4- 0.3 


0.42 


0.29 


0.29 


0.71 


0.29 


0.3- 1.0 


0.67 


0.10 


0.23 


0.77 


0.23 



Table 3: Same as Table 1, except in this case the helicity fractions are given as a 
function of cos^w- a-nd combining the y/s — 183 GeV and y/s — 189 GeV data and 
Monte Carlo. 
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Figure 1: The cos ^w- distribution for (a) W+W — > evqq and (b) W+W — > /xi/qq events from 
the a/s — 189 GeV data (points) and the KORALW Monte Carlo prediction (histogram). 
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Figure 2: Efficiency- and background-corrected cos^* distributions for (a) leptonic W decays 
and (b) for fiadronic W decays at a/s — 189 GeV. Tlie fit results for the different W fielicity 
liypotlieses are also shown. 
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Figure 3: Corrected cos 9* distribution from leptonic W decays for (a) enriched and (b) depleted 
transverse W polarisation regions together with the fit results. For (a) the required ^w- must 
satisfy 0.3 < cos^w- < 0.9, while for (b) it has to be cos^w- < 0.3 or 0.9 < cos^w 
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Figure 4: Same as Figure 3 except that in this case | cos 6* \ is shown for hadronic W decays. 
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Figure 5: W helicity fractions /o and /_ for three different bins of cos^w- in the combined 
data sample and in the KORALW Monte Carlo. 
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